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Abstract: Among the different methods to derive particle creation, finding the quantum stress 
tensor expectation value gives a covariant quantity which can be used for examining the back- 
reaction issue. However this tensor also includes vacuum polarization in a way that depends on 
the vacuum chosen. Here we review different aspects of particle creation by looking at energy 
conservation and at the quantum stress tensor. It will be shown that in the case of general spherically 
symmetric black holes that have a dynamical horizon, as occurs in a cosmological context, one cannot 
have pair creation on the horizon because this violates energy conservation. This confirms the results 
obtained in other ways in a previous paper m- Looking at the expectation value of the quantum 
stress tensor with three different definitions of the vacuum state, we study the nature of particle 
creation and vacuum polarization in black hole and cosmological models, and the associated stress 
energy tensors. We show that the thermal temperature that is calculated from the particle flux 
given by the quantum stress tensor is compatible with the temperature determined by the affine 
null parameter approach. Finally, it will be shown that in the spherically symmetric dynamic case, 
we can neglect the backscattering term and only consider the s-waves term near the future apparent 
horizon. 
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I. INTRODUCTION 

Studying Hawking radiation for cosmological black holes, which are embedded in a cosmological 
background, is not trivial [1]. One of the main problems is that due to the dynamical nature of the 
infalling matter in this case, we have to consider dynamical black hole horizons in the real expanding 
universe [2], which have different properties than stationary event horizons. 

Although about 40 years have passed since Hawking’s paper showed that Schwarzschild black holes 
emit thermal radiation, some aspects of this process are not understood totally. Hawking used quantum 
field theory in a curved spacetime to calculate the expectation value of the particle number operator 
[3] . However the expectation value of this operator does not expresses pure particle creation in a curved 
space time. We will discuss how quantum vacuum polarization may also be included in this quantity. 
Usually in the case that the spacetime is asymptotically Minkowskian in both the far past and future 
we can use the Bogoliubov transformation relating the future and past vacua to read off the pure 
particle creation in the spacetime: (0|A^fc|0) = (0|a^afc|0) = However the real universe is not 

asymptotically flat in either the future or the past. 

Hawking’s quantum field theory approach to black hole radiation, which applies to late time stationary 
black holes, is not a suitable method for calculating the Hawking temperature in the case of a fully 
dynamical black hole, where one has to solve the field equations in a changing background. There are 
alternative approaches allowing one to calculate Hawking radiation in a Schwarzschild spacetime without 
using the field equations, such as finding the related vacuum via the tunnelling method [1]. One can 
extend these methods to calculating quantum fields in the dynamical case Ell], giving us the ability 
to determine Hawking radiation for dynamical black holes in a cosmological context [Hl[5S]. However 
differentiating the different quantum effects, like vacuum polarization and particle creation, needs careful 
study of the expectation value for the quantum stress tensor {ip\T^^\'ip). 

It is known that the particle concept in quantum field theory is a global concept. The particle 
modes are defined on the whole spacetime, so that particular observers specify them by a field mode 
decomposition. There is also the number operator describing the response of a particle detector, which 
depends on an observer’s past history. To get a more usable definition of the particle state, one needs 
to construct locally-defined quantities. The best candidate for studying particle creation and other 
quantum effects locally is the expectation stress tensor value [2], which assumes a particular 

value at each point x of spacetime. This stress-tensor is objective in the sense that for a fixed state 
|^/>), the results of different measuring devices can be related in the familiar fashion by the usual tensor 
transformation rules. For example if = 0 for one observer, it will vanish for all observers. 

However this term does not express pure particle creation, it includes also vacuum polarization; and the 
way it does so depends on the choice of vacuum. 

There are several effects in quantum field theory one might take into account in a curved spacetime: 
zero point energy, the static Casimir effect, the dynamical Casimir effect (which is like a moving mirror), 
and the Schwinger effect. Particle creation can be attributed to the Hawking effect, the dynamical 
Casimir effect, and the Schwinger effect. Basically, Hawking particle creation is thermal radiation 
due to the black hole apparent horizon, which is observed at a large distance (or alternatively is due 
to an adiabatic condition between affine coordinates on past and future null infinity |24j 1: and the 
apparent horizon is the event horizon in the static case, so then the radiation is associated with the 
event horizon. On the other hand, the dynamical Casimir effect is the production of particles and 
energy from an accelerated moving mirror, so it is due to dynamical boundary conditions (see |39j for 
a comprehensive review). Finally, the Schwinger effect is a non-perturbative QED phenomenon. It 
is the spontaneous production of and e~ pairs in the presence of strong (usually constant) electric fields. 

To include all particle creation effects in the dynamical case of black hole collapse, we first need to 
solve the wave equation with suitable boundary condition in the presence of any electromagnetic field; 
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second, find an appropriate vacuum corresponding to the collapsing model (like the Unruh vacuum); 
and last, read off the general particle creation from the number density operator {tp\N = aj.afej'i/') and 
expectation stress tensor . Considering the Schwinger and dynamical Casimir particle creation 

effects are beyond the scope of this paper. Here we just discuss quantum effects due to the curved 
classical spacetime background. 

The basic point of this work is that it gives a compact review of different aspects of particle 
creation, particularly in terms of deriving the quantum stress tensor, in ways which can be applied to 
non-stationary metrics as well as the stationary case. This is what is needed for backreaction studies in 
the dynamic black hole case. 

This paper is organized as follow: We will show the relation between energy conservation for the 
particle creation mechanism and the type of the horizon in section II. In Section III, we discuss the 
particle creation concept via expectation values of stress tensor. Section IV considers particle creation 
concepts in cosmological models. In Section V we discuss different aspects of these quantum effect for 
compact stars and black hole models. We then conclude in section VI. 


II. ENERGY CONSERVATION FOR VACUUM CREATION AND ANNIHILATION 

As is known, if the total energy change for particle production does not vanish, this process is 
forbidden by energy conservation. We consider here necessary conditions so that energy conservation is 
obeyed for virtual particle pair creation and annihilation. This puts constraints on particle creation in 
dynamical black hole spacetimes. Although semi-classical physics breaks the energy conservation law up 
to the fluctuations allowed by Heisenbergs uncertainty principle AtAE = h, for long-lived real particles 
we cannot neglect energy conservation violation, so must demand AE 

It can be shown that if there is a Killing horizon, one does not have to generate energy for particle 
pair creation; thus there is no problem in the static case. But in the dynamic context, for example 
when matter and radiation fall into the black hole, there is no Killing horizon, instead there is an 
apparent horizon that is spacelike when such infall is significant [5^. We show here that this makes pair 
production inconsistent with energy conservation, and so Hawking radiation is then prohibited. However 
at late times it becomes an isolated horizon; then, one can again have particle pair creation without 
violating energy conservation. This confirms the results of pS] . 

The stationary case: Consider the possibility of particle creation by a stationary gravitational field. 
The energy of a particle in such a field is E = where is the four-momentum of the particle, and 

is the Killing vector field. The energy A of a particle is always positive outside the black hole horizon, 
where the Killing vector is timelike. The Killing vector is spacelike inside the Killing horizon = 0, 
and the energy is negative there. Therefore, this allows particle pair creation just around the Killing 
horizon. On the other hand, we know that a Killing horizon in a stationary spacetime is necessarily an 
event horizon (see Hawking and Ellis m. proposition 9.3.6). Hence, one can expect particle creation in 
a stationary spacetime which contains a black hole m- 

The dynamical case Now the question is how does energy conservation work for pair creation around 
a dynamical black hole. It is not possible to define a conserved quantity E for an evolving space time, 
as it does not have a Killing vector field. But locally for a vector field U^{x) that is timelike outside the 
black hole, we can define the particle’s energy relative to an observer with 4-velocity a.s E = —p^U^. 
For instance, in a general spherical symmetry spacetime, one can define the energy E = —p^K^ of the 
particle relative to the Kodama vector [3, which in the stationary case becomes the same as the 
Killing vector. 

Consider pair creation around the OMOTS (Outer Marginally Trapped 3-Surface) when there is 
infalling radiation or matter, so this surface is spacelike P]. Since the particle - antiparticle pair will 
be located inside the horizon (point B in the Fig. 0 ) and will each have energy defined relative to the 
same timelike vector, the energy E for both particles will have the same sign. As a result, any long-lived 
particle creation violates the energy conservation around the OMOTS. Hence Hawking radiation will 
not occur there pS] . 
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The isolated horizon case When the matter flux becomes very small, the black hole apparent 
horizon becomes an isolated horizon [19j where its tangent vector becomes a null vector and it be¬ 
comes a null surface. Note that a Killing horizon assumes the existence of a Killing vector field in some 
neighborhood of the horizon, but the isolated horizon is defined only in terms of the intrinsic and ex¬ 
trinsic geometry of the horizon itself, where is a Killing vector for the intrinsic geometry of horizon 
[19j . To consider energy conservation for pair creation around the isolated horizon, we limit ourselves to 
the general spherically symmetric dynamic black hole. Consider that the collapsing fluid within a com¬ 
pact spherically symmetric spacetime region will be described by the following metric in the comoving 
coordinates 

ds^ = + R{t, rfdn^. (1) 

The Kodama vector for this metric is 

=e-^''+^\R!,-R,G,0). ( 2 ) 

Using the Einstein equations one can show that the norm of the Kodama vector is 

= (3) 

where M is the Misner-Sharp mass for the spherically symmetric model [30) . This vector is spacelike 
inside the apparent horizon and timelike outside it. Using the energy definition E = —p^K^ , one can 
show that the energy of the particle outside the isolated horizon m is positive and the energy of the 
particle inside is negative (point A in the Fig. 0). Hence, one can have pair particle creation for isolated 
horizons without violating energy conservation. While we have shown this using the Kodama vector for 
reference, that is for convenience and is not essential; some other timelike vector field could be used for 
general spacetime. 



FIG. 1: Particle creation and annihilation around the apparent horizon for the case if an influx that decays away 
at late times. If a pair particle create at point B, they will have the same sign energy and the cannot live long 
because they violate the energy condition. But when they are created on the isolated horizon at point A the two 
particle have different sign of energy which keep the energy conservation. Note that we did not include radiation 
back reaction in this Figure, which is adapted from arXiv:gr-qc/0308033. 


III. DISTINGUISHING PARTIGLE CREATION FROM OTHER QUANTUM EFFECTS 

This section is devoted to studying the discrimination between particle creation and vacuum polariza¬ 
tion by choosing suitable coordinates, which cover different parts of the global spacetime. 

The excellent book by Frolov and Novikov m identifies (Chapter 10) three relevant quantum effects 
when discussing black holes: namely (i) particle creation, (ii) vacuum polarization, and (iii) quantum 
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fluctuations of the metric. As long as we consider a spacetime with curvature smaller than the Planck 
curvature 7^ = ^ ~ 4 * 10®®we can neglect quantum fluctuations of the metric, so we need only 
take the other two into account. 

As is known [II], quantum field theory constructed on a curved space time experiences gravitationally 
induced vacuum polarization, which is a deformation of the vacuum fluctuations by the external 
gravitational field at a given time (similar to the Schwinger deformation by an electric field). Hence, it 
is described by terms that depend only on the local curvature characterizing the gravitational field at a 
given location. This effect typically induces a nonzero vacuum expectation value for the stress energy 
tensor; it is a local effect, in contrast to the non-local nature of particle production. However the vacuum 
is not a local quantity, it must be determined from global properties; the polarization is local once the 
vacuum has been chosen. Thus because there are alternative dehnitions of the vacuum, induced vacuum 
polarization depends on the vacuum, and there is no unique way to distinguish the particle creation and 
vacuum polarization contributions to the total energy-momentum tensor. In this section we discuss how 
we can sensibly separate quantum particle creation, represented in from vacuum polarization, also 

represented in . 

Vacuum polarization and particle creation At a first glance, it seems that the expectation value 
for the quantum stress is a local quantity and hence we can straightforwardly describe particle creation 
by the local quantum stress tensor. But finding coordinates that either cover all space time, or that 
cover only the observer’s region, forces us to know global properties of the space time (even a particle 
detector that records the existence of the particle by transition to an excited state has to know about 
this interaction with the particle along the entire world line, which is non-local |35|1. For instance, if one 
has a coordinate that covers all spacetime in a case when the observer cannot see the black hole interior, 
the corresponding vacuum defined in terms of those coordinates will have a mixed state relative to the 
observer. 

If observers are sitting outside the black hole, coordinates that only cover the exterior region of the 
black hole, like tortoise lightcone coordinates (u, v) in the Schwarzschild metric, are suitable coordinate 
for stress energy examination. Here are some points that help us to specify the vacuum polarization and 
particle production: 

• If we choose a coordinate (vacuum) that only covers outside the horizon, the quantum expectation 
value of the stress tensor in this coordinate system only describes vacuum polarization. 

• One test of the vacuum polarization term is that it must be zero in flat regions of space time. 

• Another test for the vacuum polarization term is that it must be zero in local inertial coordinates. 

• For the case of particle creation, there is at least one non-zero component of the quantum stress 
tensor for all observer in the flat region of the spacetime, provided our vacuum contains the black 
hole interior region and not only outside the black hole. 

• An important test for having particle creation is that there is a non-zero flux, at the horizon. 

The two-dimensional general case We write a 2-dimensional general metric as 

ds^ = —d^~^ d^~ = —e^^dx~^ dx~ (4) 

where are local null coordinates associated with inertial coordinates, while x^ are general null coordi¬ 
nates. Notionally this is a 2-dimensional section of a 4-dimensional spherical spacetime. Here our analysis 
is for a massless scalar field, which gives the main features of the radiation (in the static case, 81 percent 
of the radiation is massless [42]). The general expression for the stress tensor in the 2-dimensional black 
hole case can be written as (see equation (5.23) in [H]): 

{^\T±±{x)\'iIj) = (V'l : T±±{x) : \^) + (5) 

where the second term is the Schwarzian derivative 

d?x^ .dx'^ 3, d^x^ ,d5:^,2 

d{x^Y dx^ 2 d{x^Y *^2;=*= 


{x^,x^} 


( 6 ) 
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A direct calculation shows this term is 

If we go to the case that a;^ coordinate cover outside the black hole (they are associated with the observer’s 
location) and the |'0) vacuum covers all space time including the black hole region, the term 2 !^ 
describes vacuum polarization, and the term ('0| : T±±(x) : |'0) gives the real particle creation part. If 
|'0) is expressed in terms of the coordinate (a;“'",a:“), the real particle creation flux relative to a observer 
with time coordinate t = ^ is equal to 


h / dx'^ 3 d^x'^ dx~^ 2 

247r \d(a::+)3 dx+ 2 d(x+y dx+ 


d^x 

d{x~)^ 


,dx 

dx~ 


3 , d?‘x 
Td(^ 



( 8 ) 


These equations are general. Their interpretation depends on the vacuum chosen. 

The four dimensional spherical case Similarly for a 4-dimensional spherically symmetric black 
hole we have (see equation (5.157) in [12] ): 

(iA|T±±(a;)|V’) = ( 1 AI : T±±{x) : IV') - ■^{d±pd±p - d\p) + ^{d±pd±(i> + p(d±4>f) (9) 

where the two-sphere radius is i? = e~^ {(p is the dilaton field) and 

ds^ = —e^^dx^dx~ + ( 10 ) 

Similar to the 2-dimensional case if the coordinates {x'^,x~) only cover the exterior of the black hole 
and the \tjj) vacuum covers all space time including the black hole region, the term 

- ^{d±pd±p- dip) + ^id±pd±p + p{d±(l))^) ( 11 ) 

is the vacuum polarization and {'ipl '■ T±±{x) : \tjj) is the pure particle creation part. Therefore the particle 
flux relative to an observer with time coordinate t = ^ is described by 

T++{x) : \Pj) - (V-| : T__(x) : |V-). (12) 

Covariantly, the created particle energy density and flux [46] are 

Pp = u^u''{tp\ : : Itp), T = : Tp^ : |'(/') (13) 

respectively where is the 4-velocity of the observer and n’' is a unit vector normal to u^. On choosing 
the 4-velocity of an observer (detector), we can find the related radiation flux and find the effective 
temperature for this observer |41|1. 

The Unruh effect As regards the Unruh effect, the normal ordered form of the quantum stress tensor, 
: Tuu ■, for an accelerated observer in the Minkowski vacuum, gives the energy that the accelerated 
detector absorbs and then makes a transition to the excited state: it is {M\ : T„„ : \M) = Here 
the tortoise coordinates (u,v) are the same as the (—,+) coordinates above. This is neither particle 
creation nor vacuum polarization, because it is in Minkowski spacetime. It is just an excited state due 
to acceleration (a ’’fictitious particle”). 

The Schwarzschild case The case of a static black hole is gievn by Schwarzschild space time and 
its vacuums. There are two well known null coordinates for the Schwarzschild metric. The first is the 
tortoise lightcone coordinates 

u = t — r*, v = t + r* (14) 

where r* = r — 2m + 2mln(^ — 1). These coordinates for r > 2m vary on the range —00 < u,v < -|-oo. 
The second coordinates are Kruskal-Szekeres null coordinates {U, V) defined by 

U V 

U = —4m exp(—-—), V = 4m exp(-—) 

4m 4m 


( 15 ) 
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There are three vacuums for the Schwarzschild metric, which are defined as below [22j : 

(i) the Boulware vacuum \B), defined by requiring normal modes to be positive frequency with respect 
to the Killing vector with respect to which the exterior region is static. This vacuum is defined only 
in the exterior region. The expectation value of stress tensor in this vacuum only describes vacuum 
polarization. 


(ii) the Unruh vacuum |t/), defined by taking modes that are incoming from to be positive 
frequency with respect to while those that emanate from the past horizon are taken to be positive 
frequency with respect to U, the canonical affine parameter on the past horizon. This vacuum is defined 
on the exterior region together with the interior region associated with particles infalling from that 
exterior region. The expectation value of the stress tensor for this vacuum contains both particle creation 
and vacuum polarization parts, and has a non-zero flux |21j . 


(iii) the Hartle-Hawking vacuum |i7), defined by taking incoming modes to be positive frequency with 
respect to V, the canonical affine parameter on the future horizon, and outgoing modes to be positive 
frequency with respect to U. This covers the entire maximally symmetric vacuum Schwarzschild solution. 
The expectation value of the stress tensor in this vacuum contains both particle creation and vacuum 
polarization parts. This has zero flux (12), showing this is thermal equilibrium of a black hole with its 
emitted radiation. 
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FIG. 2: The expectation value of the energy momentum tensor for different vacuums in 2-dimensional 

Schwarzschild spacetime. The left number in each column was calculated at the horizon and the right number 
is at the infinity. The last column represents the radiation flux. By equation 0. vacuum polarization is given by 
subtracting the third column from the first, and the second column from the fourth (there is another component 
not shown here). 


A summary of the expectation values of the energy momentum tensor for 2-dimensional Schwarzschild 
spacetime is shown in the Table 

Conceptually, the radiation flux appears for an observer that cannot see parts of the spacetime. Since 
the Hartle-Hawking vacuum and Unruh vacuum cover inside the horizon, which can not be seen by 
external observers, the related states for these vacuums are mixed states. 


IV. PARTICLE CREATION IN COSMOLOGICAL MODELS 

In this section we investigate particle creation in FLRW metrics using the quantum stress tensor. 
In particular we discuss de Sitter spacetime in stationary and cosmological coordinates and define the 



























related vacua for these cases. Other aspects of radiation production like the affine parameter approach 
and the relation to the inflationary epoch are discussed. 

Particle creation in cosmological models is substantially different from that in black hole models. The 
key difference is that in the cosmological case, if event horizons exist for an observer, which will be the 
case for de Sitter spacetime and for expanding universes with a cosmological constant, every observer 
can at very late times receive thermal radiation associated with their cosmological event horizon (note 
that different observers will have different event horizons). Hence, thermal radiation will in principle 
appear for observers everywhere. However for the black hole case, there is a unique event horizon for all 
external observers and only a distant observer outside the black hole will see such radiation at late times. 
Therefore, we have well-defined observers in the black hole case for whom the thermal temperature and 
backreaction can be calculated, but in the case of cosmological models, defining the temperature and 
back reaction of the radiation has ambiguities that must be solved. 

One should note here that there are various calculations of the vacuum energy of zero point quantum 
fluctuations, which in principle is one of the main potential candidates for cosmological dark energy 
[Min]- This zero-point energy of a quantum field diverges, and so is calculated up to a cut-off scale 
in the case of both flat and curved spacetimes, and turns out to be hugely larger than the observed 
cosmological constant, which is a major problem for cosmology if this field gravitates (which need not 
necessarily be the case [H])- This however is an essentially local effect that can occur in cosmology as 
well as in the black hole case, but is not directly related to either particle creation or vacuum polarization 
effects, which are confined to curved spacetime backgrounds. We will not consider this further here. 


The standard model in cosmology is based on a homogeneous and isotropic expanding universe. The 
metric which best describes the early accelerating expansion (inflation era) and late time acceleration 
phase (dark energy dominated) is the de Sitter spacetime in non-stationary coordinates, given by the 
metric (35). These coordinates only cover half of the de Sitter manifold — -\-Z\-\- Z'^ = A > 0. 

However because it is a spacetime of maximal symmetry, one can also use static coordinates with the 
metric given by (16). These also only cover part of the de Sitter manifold. The matter dominated 
and radiation dominated stages of the universe after inflation are well described by a flat Friedmann- 
Lemaitre-Robertson-Walker (FLRW) model. In contrast to the de Sitter spacetime, this does not have a 
cosmological event horizon. 


A full cosmological model will have three epochs: an initial inflationary era with metric close to de 
Sitter, followed by a matter and radiation dominated era, and finally a late time dark energy dominated 
era that is again like de Sitter. In this section, we discuss particle creation in these various metrics. 
But one should note the following: the event horizon in cosmology only refers to the far future of the 
universe. If creation occurs relative to that horizon, it has nothing to do with anything we can observe 
at the present day. 


A. Particle creation in FLRW spacetime 


The first prominent work to calculate particle creation in a FLRW metric is Parker’s work in 1969 
[31], which calculated the probability of finding a particle at the time t according to the Bogolubov 
coefficients. Being a conformally flat spacetime, one has the ability to specify a well-defined confor¬ 
mal vacuum and determine the mean particle creation rate in conformal time. As discussed above, 
such non-zero particle creation, which is calculated by comparing the different vacuums according 
to the Bogolubov coefficients, cannot be considered as pure particle creation: it can be fictitious or 
quasi-particles. To study the real particle creation rate, one needs to determine that rate ‘from (^|T^y |'0). 


Using equation (A4), it can be shown that the radiation flux (12) for the minimally conformally 
coupled massless sca ler f ield is zero for 2-dimensional FLRW models. For 4-dimensional FRLW metrics, 
using the equation (B9| to calculate {Ttx)-, one can show that the radiating particle flux is zero for 
conformally coupled massless scalar fields. The non-zero component from (A4| and (B9) can be 
interpreted as vacuum polarization. 


One can define a cosmological apparent horizon as follows. iJ is a hypersurface in a 4-dimensional 
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spacetime that is foliated by 2-surfaces such that \h= 0, 0(f) \h^ 0, and £e9(n) \h^ 0. This 
separates regions where either future or past null geodesics going either both inwards or both outwards 
converge, from where this is not the case. An apparent horizon is called outer if £g9(^n) \h< Oj inner 
if £i9(^n) \h> 0, past if 0(f) \h< 0, and future if 0(f) |fr> 0. They are trapping horizons, that is, they 
are related to spacetime regions where particles are causally trapped in the sense that they cannot 
reach infinity, if they are future directed, but not if they are past directed. Past directed apparent 
horizons occur in all realistic cosmologies and are related to the existence of the singularity at the start 
of the universe nnj. Future directed apparent trapping horizons occur in cosmology only in universes 
that recollapse to a second singularity in the future, when they are related to the existence of that 
singularity. The most relevant case in the context of an ever-accelerating expanding cosmology is the 
inner future apparent horizon (IFAH) which can occur if the cosmological constant is positive. The 
physical intuition for this horizon is that it is the boundary of the region so that to the future, the 
cosmological expansion is so strong that even the future directed ingoing null geodesic cannot con¬ 
verge. For de Sitter spacetime this cosmological horizon coincides with the de Sitter event horizon. It is 
not a trapping horizon because it is future directed; however it can be associated with Hawking radiation. 

Recently some workers have applied the particle tunneling method to show that there is particle thermal 
radiation from the FLRW apparent horizon (cosmological horizon) in the dynamic case [27] . They have 
neglected a basic point. Their calculation is based in the WKB approximation assumed for the light 
coming from the apparent horizon to the observer. As discussed in [5SJI15], the WKB approximation 
cannot hold for light passing the FLRW apparent horizon, except in the case of de Sitter spacetime. As 
a result, there is no particle creation in the FLRW metric in the matter and radiation dominated eras. 


B. Particle creation in de Sitter spacetime 

Here, we consider particle creation in the stationary and cosmological versions of de Sitter spacetime. 


1. Particle creation in the stationary de Sitter spacetime 


Historically, by using the path integral method Gibbons and Hawking showed that any observer 
will see thermal radiation in de Sitter spacetime. Let us start our calculation from de Sitter spacetime 
in stationary coordinates: 


ds^ = -(1 - ^R^)dT^ + {1- ^R^y^dR^ + R'^dn^. 

o o 


(16) 


This metric has a time like Killing vector field d/dT, and R = ^ is a comoving observer’s event 

horizon with surface gravity Kc = These coordinates cover only part of the de Sitter hyperboloid 

[To] : they cover the triangle on the left in Fig.([^ that stretches from to *+. 

The null coordinates Defining the new variable r*(r) by dr* = 


2g4 


In 


l+yfi?' 


(17) 


(18) 


These coordinates are defined for 0 < i? < y vary over the range —oo < u,v < -|-oo. The second 

coordinate definition is Kruskal-Szekeres-like null coordinates {U, V) given bv[47| 


Similar to the Schwarzschild coordinate we define tortoise lightcone coordinates (u, v) by 

u = t — r*, V = t + r*. 


^ = \/T®^P(^cu), V =-\ -eyip{-Kcv). 


(19) 
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In this coordinate system, the metric becomes 

ds^ = -^- dU dV + R^dn^. (20) 

As shown in Fig.0, light coming from past null infinity, , and passing near the event horizon is 

stretched infinitely. The different with the black hole case is that here we must consider ingoing null 
geodesics instead of outgoing ones. The observer coordinates (u, v) cannot see outside the event horizon 
while the (U, V) coordinates cover both outside and inside the horizon. 


The various vacua Like the Schwarzschild case, the stationary de Sitter spacetime has three vacua. 

(i) We define the \B) vacuum state (which is like the Boulware vacuum) for the wave solution in terms 
of the mode functions {u,v). The corresponding modes are 


1 


1 


•s/dTl 




( 21 ) 


It is not suitable as a vacuum for the stationary de Sitter spacetime with particle production because it 
covers only the inside of the horizon. 

(ii) We define the \H) vacuum state (which is like the Hawking-Hartle vacuum [37]) for the wave 
solution in terms of the mode functions (U, V). The corresponding modes are 


^—iwU 


—iwV 


•s/dyi 


•v/Di 


( 22 ) 


It is not suitable as a vacuum for the stationary de Sitter spacetime because it covers the inside and 
outside of the horizon in a time symmetric way. It can be shown that 


(H\ : T,,(x) : \H) = (H\ : T„„(a;) : \H) = 


hA 

1447r 


(23) 


at the center, and hence the matter creation flux is zero. It describes a thermal equilibrium universe 
with temperature T = ^. 

(iii) We define the \U) vacuum state (which is like Unruh vacuum ) for the wave solution in terms of 
the mode functions (u,V). This mode function only covers the upper triangle in the Fig. @ , and it can 
be a suitable vacuum for particle production in the stationary de Sitter spacetime because it covers both 
inside and outside the horizon in a time asymmetric way. The corresponding modes are 


^—iwV 


^/T 


TTW 




TTW 


(24) 


As we attribute the \U) vacuum to the {u,V) coordinates, we can calculate the matter creation flux 
from equation ®: 




where the thermal temperature in this case is 


T = 


HKc 


(25) 


(26) 


From the energy conservation point of view (Section II), a Killing horizon exists, the Kodama vector 
changes sign on the event horizon, and hence energy conservation is valid for the created pair. 


Affine parameter approach We can also derive this temperature from the affine null parameter 
approach [24]. The metric in this case is stationary and the adiabatic condition is satisfied for null 
geodesics which come from near the event horizon. Assume that we consider the upper triangle in Fig. (§, 
which is covered by the de Sitter space coordinates (V, u). It can be shown that the v parameter is the 
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FIG. 3: Particle creation in the stationary de Sitter spacetime 


affine null parameter at future null infinity where the observer receives the radiation, and V coordinate 
is the past null infinity affine parameter. Therefore, we can write the equation 


V = p(v) 



exp(—Kcv) 


(27) 


for the null geodesic which comes from past null infinity and arrives at the observer’s point in future null 
infinity. The adiabatic condition will be held for the light passing near the de Sitter space time horizon 
pS] and it can be shown that the spectrum of the radiation is Planckian (23]. Assuming the adiabatic 
approximation for the null geodesic near the event horizon, we get 


1 p(y) 
2tt p(v) 


(28) 


Using the equation (19) we get the temperature T = ^. This temperature is consistent with the 
temperature (26) from the expectation value of quantum stress tensor (25). 


Note however a key feature this derivation makes clear: the observer only experiences particle creation 
at the end of her history where relation (27) diverges, and hence will not experience it at any finite time 
before then. Hence static de Sitter spacetime regions are not filled with black body radiation everywhere: 
it comes in to being for an observer only at very late times. 


Vacuum Polarisation The pure vacuum polarization effect for the 2-dimensional de Sitter spacetime 
can be calculated in tortoise coordinate (u, v) by the terms 

(H|r„„|H), (H|T„„|H), (H|r™|H), (29) 

where the \B) state is like the Boulware state for the tortoise coordinate. The de Sitter metric in the 
tortoise coordinates is 


ds^ = -4 




(1 


du dv. 


(30) 


Here comparing with the equation (4) we get p = In 2 


e2 V 3 


(u-„) 


l+eV 3 




According to equation (7), the 
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different components of the vacuum polarization are, 

{B\T^^\B) = -^{dupdup-dlp) 

“ 3 ^2 " 1 + ev^ («-«) ^ ^ + j “ "144^ 

(BIT,, IB) = 

/i /a 1 eA(“-’') A \_ M 

127r 3 2 I ^\/^{u-v) 3 I ^y^iu-v) j 1447r 

A esAl"-*') \ _ HA 

3^i + eV¥(“-r ; 


(B|T,„|B) = -A(a^pa,p-Cp) 

127r y 3 2 X _|_ gV^(“-*') 


(31) 


(32) 


(33) 


These equations show that the vacuum polarization flux relative to a observer with time coordinate 
t = is zero. The vacuum polarization is homogeneous for the de Sitter spacetime, and it reduces to 
zero in the Minkowski limit A —> 0, which show these results are due to gravitational effects in quantum 
field theory. The energy density of the vacuum polarization for the observer who is sitting in the center 
r = 0<t4>u — u = 0is zero: 


Pvp — (B|r„,|B) + (B|T,,|B) + 2(B|T„„|B) — 0. 


(34) 


The total effects A summary of the expectation values of the energy momentum tensor for 2- 
dimensional static de Sitter spacetime is given in Fig.(^. Here we see that particle creation flux occurs 
only for the Unruh-like vacuum and the expectation value for the Boulware vacuum state represents pure 
vacuum polarization. 


''^vacuiir^^ 

(7-uu) 

<rw) 

{:7-uu;> 

(; 7-w:) 

f=(;rw:)-{;ruu:) 

Boulware (u,v) 

AA 

" 14471 

AA 

' 144 77 

0 

0 

0 

Hawking-Hartle 

0 

0 

ftA 

ftA 

0 

(U,V) 



144 71 

144 71 


Unruh-like (u.V) 

AA 

0 

0 

ftA 

ftA 





144 71 

144 71 


FIG. 4: The expectation value of the energy momentum tensor in different vacuum for 2-dimensional de Sitter 
spacetime observed at the center. The meaning is as in the case of Fig. 
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2. Particle creation in the cosmological de Sitter spacetime 


The cosmological de Sitter spacetime is described by the metric 

ds^ = -dt^ + e^^^idr'^ + r'^dVi^). (35) 

These coordinates only cover half the de Sitter manifold —Zq + = I-, A > 0; consequently 

as shown in Fig.([^, this spacetime covers the upper half of the total de Sitter manifold. Like the case of 

stationary coordinates, this spacetime has an event horizon at r = and the {t,r) coordinates 

cover both inside and outside the horizon. Using the coordinate transformation, 


R = 


3-r, 


T = t 


3 A 

4- dR = t - — In (1 - -i?2) 

1 - 4i?2 2A ^ 3 ^ 


(36) 


we can relate stationary de Sitter coordinates (T, R) with metric (161 to the cosmological coordinates 
(t, r) with metric (35). 


This coordinate transformation changes the particle creation effects discussed above in the static case 
for an observer who is located at r = 0 in the expanding case. Particle creation for the cosmological de 
Sitter spacetime can be seen by other methods as well PEg. Again the effect only occurs at late times 
when the affine parameter relation diverges. It does not manifest for an observer at a hnite time in her 
history (because no matter what that hnite time is, there is still an inhnite time to run before the event 
horizon is reached). 


Vacuum Polarization Using equation (A4|, the expectation value of the quantum stress tensor for 
a general FLRW metric is 


{Tun) 


^ h{ad - a?) 
' 487r 


—H{ad) 

487r 


(37) 


This shows that the FLRW model has no hux of Hawking radiation and no vacuum polarization. Thus 
in the case of de Sitter spacetime we get 


(Tun) = (T™) = 0, (T„„) = (38) 

We can calculate the energy density of vacuum polarization(VP) for the comoving observer = (1,0) 
inside the horizon as: 


Pvp = -^i{x'^\Tuu\x'^) + (a:=^|r^„|a;=^) + 2(a:=^|r„„|a;=^)) = 


—hA 
72 tt ■ 


(39) 


Like in the case of stationary coordinates, the vacuum polarization in a cosmological de Sitter spacetime 
is homogenous. The comoving cosmological observer will see a constant vacuum polarization energy 
density which is smaller in magnitude than the cosmic fluid by order h . 


The Kodama vector ([^ changes sign on the event horizon, which is a null surface, so particles and anti¬ 
particles have different signs for their energy and hence energy conservation is obeyed for any particles 
created there; hence particle creation is possible. One might think that such particles created at the event 


horizon are diluted away by cosmological infinite redshift, 1 -I- Zc 



but this argument is not 


true. This is like the black hole case where the created particles near the horizon get infinitely redshifted 
when it reaches future infinity. 


3. The expanding de Sitter era in standard cosmology 


The standard model of cosmology uses the expanding de Sitter metric (35) for describing both the 
early inflationary era and the late dark energy dominated era, but with different values of the constant 
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A in these two phases, which are separated by a matter and radiation dominated era. Since there is 
both a start time and an end time for the inflationary era, this phase of the universe’s history is only 
covered by the small part of the de Sitter spacetime which is located between times ti and t/ in Fig.([^. 
The late dark energy dominated phase starts at a time toE and (assuming it is a cosmological constant) 
lasts forever, and so has an associated event horizon as t —>■ oo. 

Late time particle creation Observers have an event horizon when vacuum energy dominates at 
late times (so future infinity is spacelike) so one might expect Hawking radiation to be associated with 
this horizon at very late times. Note that the event horizon for late time accelerating phase is the same 
as the cosmological apparent horizon. Overall, the initial state (the vacuum choice at the start) and 
state at the event horizon of the universe need knowledge of the pre-inflation era and final fate of our 
universe respectively. If we have a finite time for an exponentially accelerating era, it is hard finding the 
particle creation term from the quantum stress tensor when we do not know the initial state. We can 
however use the tunneling method [6] to calculate the flux of radiation at very late times. However since 
any particles created near the cosmological event horizon would arrive at the center at very late times, 
in practice a comoving cosmological observer cannot see any associated created particles at the present 
time in an expanding universe. They will arrive at her world line only in the very far future. Thus 
particle creation in the late dark energy dominated phase will have no present day observational ef¬ 
fects in cosmology, because the event horizon for any observer only exists in the far future of the universe. 



FIG. 5: Particle creation in the cosmological de Sitter spacetime 


Inflationary particle creation Note that the de Sitter horizon for the inflationary era is not a 
cosmological event horizon, because the inflationary era comes to an end and is replaced by radiation 
dominated era. There is however a cosmological apparent horizon (an anti trapping surface), and this 
can lead to particle creation (because Hawking radiation can be associated with apparent horizons rather 
than event horizons 0 ). 

To study particle creation in the inflationary era, we have to notice these points: 

• The standard choice of the vacuum state is the Minkowski vacuum of a comoving observer in the 
far past t = —oo (when all comoving scales were far inside the Hubble horizon), which is called the 
Bunch-Davies (B-D vacuum), and is shown in Fig.([^. There are some discussions about suitability 
of choosing this vacuum |32| because, it is defined at past time infinity, which is not applicable 
to standard big-bang cosmology where inflation started a finite time ago. Changing the vacuum 
will result in a different power spectrum of the quantum fluctuations and expectation value of 
the stress tensor. For example as we saw in the last subsection, if we choose another vacuum like 
the vacuum related to the maximally extended de Sitter spacetime, there is no matter creation 
flux. Even observational parameters like the primordial power spectrum and bispectrum in the 
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CMB change on choosing another vacuum |38j . To know the correct vacuum, we have to know 
the pre-inflationary models to calculate the initial state. Since the initial state of the big bang 
cosmology is directly connected to quantum gravity problems, we cannot know the initial state at 
present. It seems that effective field theory can be a tool to estimate the effective initial state for 
pre-inflationary models. 


• The un-renormalized expectation value for the stress tensor has ultraviolet (UV) divergences in 
the B-D vacuum which cannot be renormalized by the usual counterterm. Recently, people have 
applied different renormalization methods for UV divergences in cosmological models [31] . This UV 
divergency comes from choosing the B-D vacuum, or is generated by the unphysical procedure of 
sudden matching between different cosmological eras [33] . 

• The expectation value of the quantum stress tensor in the B-D vacuum according to the observable 
coordinates describes vacuum polarization, particle creation, and vacuum energy effects (in the 
case that we consider background wave number renormalized by the cut off), see |33| . As shown in 
Fig.Q, choosing a vacuum like the B-D vacuum, which represents the initial state of the field at 
far past infinity (deep inside horizon), generates a particle creation term in the expectation value 
of the field in inflationary models. However inflation starts from an initial time, not from infinity. 


• When we change coordinates from the cosmological frame to the stationary frame, the cosmological 
redshift changes to a gravitational redshift because of that transformation. The point is that the 
attributed temperature (or wave frequency) is calculated by comparing the vacuum at future and 
past infinities, giving the temperature at the observer’s point in the future. To see this, consider the 
Kodama vector given by equation ([^. For the stationary de Sitter spacetime this is (—|Ai?, 0, 0,0) 

and for the cosmological de Sitter spacetime it is (1, — 0, 0)). One can argue that Kodama 

observers who lie inside the horizon and have 4-velocity vector K = see the divergence in 
frequency measured by such observers at the horizon u = . Therefore, such observers measure 

a thermal spectrum with temperature T = (see [8] for more discussion) which has a divergence 
at the cosmological apparent horizon. 

• An applicable method for this era is the tunnelling method. The advantage of this method is that 
we do not need to know the initial state of the pre-inflation era or the global event horizon. We just 
need the WKB approximation [28| for ingoing waves near the de Sitter horizon (the cosmological 
apparent horizon). In contrast to the black hole case where we have a apparent horizon for outgoing 
null geodesics, which is a trapping horizon, here we have a cosmological apparent horizon for ingoing 
null geodesics, which is an anti-trapping horizon (it is a time reversed trapping horizon). In black 
hole cases, outgoing null geodesics are received by the observer; here ingoing null geodesics are 
received by the observer. For the case of a de Sitter era that comes to end we have a finite burst of 
Hawking radiation flux associated with the apparent horizon which reaches the observer at very late 
times and is diluted by cosmological expansion by a factor (y^)^ (the reason for power 2 is that 
one factor is due to the cosmological Doppler redshift of the photon and second is because of 
the arrival rate of the photons is also reduced by the same factor) where z is cosmological redshift. 
The resulting Planckian spectral flux is 


p{w) 


dw r('u;) 
27r - 1 


(40) 


where F(w) is the frequency-dependent (graybody) transmission coefficient for the outgoing 
particle to reach future infinity without backscattering. The total luminosity for this radiation 
is F = p{w). As a result of the redshift factors, this particle creation term does not give a 

significant contribution to the total matter in the universe or to its expansion rate. Nor does it 
affect any observations we might make at the present time. 


But in any case, one does not need any concept of horizons or particle pair creation in order to derive 
the quantum fluctuations that lead to structure formation in an inflationary universe: see e.g. [43]. 
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The basis of inflationary model perturbations is a mechanism that generates quantum fluctuation of the 
matter field because of the non-commutativity of the field state and momentum; these become classical 
after horizon exit and then become seeds for cosmological structure formation |43j . 


V. DIFFERENT ASPECTS OF SEMI-CLASSICAL QUANTUM EFFECT FOR COMPACT 

STAR AND BLACK HOLE 

In this section, we discuss different aspects of quantum vacuum, like vacuum polarization, for compact 
stars and asymptotically flat spacetimes. For the black hole case we consider the effect of applying the 
adiabatic condition on the quantum stress tensor, looking at the spectrum of negative energy particles 
which fall into the black hole and considering the back-scattering of radiation. 


A. Quantum stress tensor for compact stars 

If the spacetime does not possess an event or apparent horizon (the case of stars and planets) then the 
above complications go away and we only have one vacuum state to deal with; the Boulware vacuum, 
which encompasses the entire space. There is no event horizon, and so there are no processes associated 
with an event horizon. 

The vacuum solution around the compact star can be described by a Schwarzschild exterior solution 
in tortoise coordinates: 


= —(1- )dudv + r^d^^ (41) 

r 

where r > 2m. In these coordinates, the vacuum polarization for the compact star is the same term 
{B\Tfj_,^\B) as discussed before [20]. One can take the vacuum solution before the star collapses, 
\in), and calculate the quantum stress tensor after the star collapses according to this prescription. 
On the other hand we can read the Bogoliubov coefficients and find the particle number operator, 
{in\Ni\in) = Both of these methods calculate the same vacuum 

polarization effect (not particle creation). Note that Fabbri et al. [33] have shown that the particle 
number operator is proportional to total quantum energy in 2-dimensional space time. 

As an example, consider a static star with spherically symmetric distribution of perfect fluid. The 
line-element can be written. 


ds^ = -e^°‘^^'>df + + r^dn^. (42) 

One can define double null coordinates (m, v) for this metric by du = dt — dr* and dv = dt + dr* where 
dr* = For these coordinates, or any other null coordinate {U{u), Id(u)), we have fy = ~fy- 

Therefore from the equations Q and 0 the vacuum polarization components give the equation 
(T„„(x)) = {Tyy(x)). Therefore, like in the Schwarzschild case, the flux of vacuum polarization for a 
comoving observer is zero. As a result, in the static spherically symmetric metric the general matter 
flux relative to the comoving observer, given by = JVp + -^pCj is equal to only the particle creation 
flux, i.e. Bvp = 0. But there is no particle creation flux as there is no event horizon for compact stars 
so Bpc = 0 = 0. 


If a central object in asymptotically flat spacetime is sufficiently non-compact, 2m <C r, then the 
spacetime is closely Minkowskian and there are no horizons.In this case the quantum vacuum polarization 
effect is proportional to the two factors h and ^ [30] [33] . The first factor shows that this effect is a small 
quantum correction. The second factor shows that such a vacuum state will be virtually indistinguishable 
from the Minkowski vacuum state. Hence, the expectation value of the (renormalized) quantum stress 
tensor will be negligible throughout the entire spacetime. In other words, this case is the weak field limit, 
and since the linear relation between two coordinates and is dominant in the equations Q and 
0, the vacuum polarization term is zero. 
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B. Black Hole Spacetimes 

Here we want to find the trace of the adiabatic condition |24] in the quantum stress tensor when black 
holes occur. Let the affine parameters on the null generators of past and future null infinity, and 
, be U and u. Null geodesics from past null infinity pass through the center and reflect off the center 
at r = 0. We express the relation between these two affine parameters provided by the null curves as 
U = p{u) and define the function 


^ P{u) 

p(u)' 


(43) 


where ' = With this definition we can show that = —k + The non-zero term in the particle 
flux equation § in the two dimensional case becomes 

(V'l : T__(x) : IV^) = y) (44) 

where U = x~ and u = x~. Now assuming the adiabatic condition |k(w*)| <C along the m* 

reference null curve completely traversing the body, we get a non-zero term for particle flux. For the 
stationary case k = so we get the standard thermal flux 


:r__(x):|V^) = 


h 


7687rM2' 


(45) 


This calculation can be extended to 4-dimensional spacetime using the equation (121 (equation (5.155) 
in [H]). As discussed in |21], the adiabatic condition (eikonal approximation) is equivalent to the state¬ 
ment that a photon emitted near the peak of the Planckian spectrum should not see a large fractional 
change in the peak energy of the spectrum over one oscillation of the electromagnetic field (that is, the 
change in spacetime geometry is adiabatic as seen by a photon near the peak of the Hawking spectrum). 
This condition on the tunneling method for black holes radiation gives a constraint on black holes that 
can emit radiation 1251. 


C. Black hole radiation from quantum stress tensor 

Let us discuss about the following important question: 

Can we recognize particle creation from the local character of the quantum stress tensor, or 
do we need some global information about the spacetime? 

Generally the final fate of gravitational collapse is a black hole or compact star. For compact stars, 
only the vacuum polarization contributes to the quantum stress. For asymptotically flat or conformally 
fiat space time like FLRW, since the affine parameters Um (we drop the in index in the formulas) and 
u on the null generators of past and future null infinity are linearly related, the particle creation flux 
part is zero from the equation ([8). But gravitational collapse that has its final fate as a black hole has 
two affine parameters exponentia ly related near the horizon. This exponential relation gives a non-zero 
term for the created matter flux. 

Collapse to a Schwarzschild black hole: Event Horizons For more clarification, consider collapse 
to a black hole so that a Schwarzschild black hole forms finally. Following the calculation in [T^ (equation 
(4.24)), if A is an affine parameter which separates outgoing null geodesic at future null infinity and ingoing 
null geodesics at past null infinity, and u = t — r*, we get 

u = 2EX — (46) 

fci 


at future infinity and 


u-Uo = -hx 


( 47 ) 
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at past infinity (ki and ^2 are constants). Therefore, we can write 

u = -‘^{U - Uo) - 4m\n{- ^ ^° ) 

k,2 kik2 


(48) 


The quantity Uq determines the last ingoing null ray that can come out of the black hole, as shown 
in Fig.( |^. Far from the horizon, the relation between two afSne parameter is linear, and according to 
relation (llj) there is no particle creation. In addition, in the zero-limit for the black hole mass we have 
zero particle creation. Only near the black hole horizon do we get the exponential relation between 
parameters and on using the equation Q we get a non-zero particle flux. Using the affine relation (48), 
we get the Bogolubov coefficients : 


= -c 


/ O / 

dsi-) 

-oo W 


1 i(w'-w^) 

2 g ' fco •' 


*=2 exp[iw4m\n{s/kik2)] 


(49) 


and 


= -c 


(50) 


ds( —) 2 e *^“''=2 expbw4TOln(—s/fciA: 2 )] 

J-oo w 

(this is an extension of equations (4.75) and (4.76) in [E]). 

In the case ^ exp{8'!Tmw)\j3ww'\‘^ we have deviation from thermal radiation at a distance from 

the horizon. All observer paths will end at as shown in Fig.(|^. Hence, all observers will see thermal 
radiation coming from near the horizon at late times. Note that all the above calculation are related to 
the case when we know the final (global) scenario to be collapse to a black hole and can use equation 
(48) for the Schwarzschild geometry. If the final fate of the collapse is a star, the linear relation between 


affine parameters gives no particle creation. We are also here assuming no infalling matter or radiation. 
Note that one can calculate the expectation value of the stress tensor in terms of the past null infinity 
affine parameters (Ui„,Un)j which are related to the future affine parameter by equation (48). This 
vacuum state is called the ”in” vacuum state, and the resulting expectation values of stress tensor are 
calculated for the Schwarzschild metric in the chapter 5 of the book [T2] . 


Collapse to a Schwarzschild black hole: Apparent Horizons For calculating particle creation 
at the point O on the apparent horizon (MOTS), we need to consider two ingoing null rays c and d and 
their affine distance at past null inhnity and at the future singularity (since these two null rays end at the 
singularity). We compare the affine distance of these two null rays before falling into the singularity with 
their afhne distance at past null inhnity. Since we are inside the event horizon the t and r coordinate 
change their character. We dehne new coordinates ^ ~ ~ ~ ^) u’ = r* — t. Then the affine 

distance relation gives 


u'= 2UA -I- 4mln(-^) (51) 

ki 

at future inhnity where E < 0 and fci, A > 0 and 

U-Uo = k2\ (52) 

at past inhnity where ^2 > 0. Therefore, we can write 

u = - Uo) -k4TOln(- ^ ^° ). (53) 

^2 kik2 

We see that similar relations will hold between the afhne parameters. With the same reasoning as in the 
case of outside the horizon, the second exponential part gives a non zero term for particle creation in 
equation 0 , which is signihcant for outgoing null rays which pass near the event horizon. The important 
point is that the — sign behind the second term in the above equation causes that the thermal radiation 
near the horizon but inside the horizon has negative energy. Calculation of the Bogoliubov coefficient 
near the horizon and inside the black hole gives 

|awP = e-8™"'|/3wP (54) 
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Hence the number of created particle at late time near and inside the horizon is 

(") = (55) 

where r(u') is backscattering factor for wave packet. 



FIG. 6: Null trajectory for the spherical black hole collapse 


Since the expectation value of the quantum stress tensor has an internal integration over the wave 
frequency, we cannot read the spectrum of the radiation from it. As is discussed in [21], the adiabatic 
condition is the necessary physical condition for having a Planckian spectrum. For a collapsing model 
that terminates in a Schwarzschild black hole, the Planckian spectrum comes from the logarithmic term 
(48) which means that only waves that pass near the event horizon have a Planckian spectrum. 


D. An example: two dimensional Schwarzschild black holes 

As is known, the tensor {T^^) has two advantages in determining particle creation effects: the first is 
that it gives a covariant definition of real particle creation , and the second is that if it becomes zero in 
one coordinate system, it remains zero in all other coordinates. To calculate the quantum stress tensor 
in the Openheimer-Snyder collapse case, consider the interior FLRW metric as 

+ r^dn^'^ , (56) 

where a{t) is the scale factor, k is the spatial curvature, r is the comoving radius, and t is the proper time 
of the comoving observer, that is, the cosmological time. The double null form of the metric is given by 

dsp =—a^dudv + , (57) 

where 

dv dx ■■= . , (58a) 

a V1 - kr^ 

u ■■= V-X, v:=r] + x- (58b) 

Assume the light which comes from past null infinity where the metric is Minkowski space time passes 
through the center and reaches a point inside the FLRW space time. Let u = x~ and v = x'^ be a null 


dsp = —dt^ + a{tY 


dr^ 


1 — 
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coordinate inside the FLRW part and u = x and v = be initial null coordinates at past null infinity 
which we denote by|a:^) and |a;^) . As shown in [12] 


{x |T±±|a; ) = {x |T±±|x ) - ,x } 


The observer’s radiation flux at the observation moment tjq is given by 

(x=‘=|T^|x=‘=) = \{x^\Tr,x\x^) = - r^^)|x=‘=) 


Substituting equation (59) into the above equation one gets 

(s±R|s±)= 


^ cPu ,du\ 

247ran V 2 du^ du du^ du J 


(59) 


(60) 


(61) 


Without giving the details of the calculation, the u which comes from the Minkowski part at past null 
infinity and passes through the center to arrive at an IMOTS (the Inner Marginally Trapped 3-Surface) 
point inside the FLRW part varies as 


U = CU + Cl 


(62) 


where c and ci are constants. As a result {x^\T^\x^) = 0, which says no quantum particle creation occurs 
from the IMOTS in 2-dimensional OS black hole. Any other collapsing models which keeps the linear 
function (62) between the two affine parameters also gives a zero term for quantum particle creation. 
Note that we have ignored the reflection of the star surface. This result is compatible with the result at 
[25] where the authors show that by the tunneling method that the IMOTS does not have any radiation. 


E. Back scattering in dynamic black holes 


In this subsection, we will show that like in the case of the Schwarzschild stationary black hole, the 
2-dimensional dynamic black hole and s-wave approximation for 4-dimensional dynamic black hole can 
describe the main features of black hole radiation. 


We write the metric in the form ds^ = —2f{u,v)du dv + R{u,v)^dn^ where u = x~ and v = x'^ only 
cover outside the horizon where the observer is. Without giving the details, the Klein-Gordon equation 
for a general wave equation becomes 




Defining the function (p = ^, we get this form for the general wave equation: 




R 


i?2 


(J) 


2i?3 


-$ = 0 . 


(63) 


(64) 


We call Vi = — ^+^2 ^ effective potential in the general spherically symmetric spacetime, 

and the wave equation reduce to this simple form 

(^+W)$ = 0. (65) 

H 

In the special s-wave case, £ = 0 and the Klein-Gordon equation is 

d+{R^d-(p) + d-{R^d+<p) = 0 (66) 


which can be written: 


d+d- 


ip- 


[d+Rd-ip + d-Rd+p>] [9+d-^p + 9_ d+ip] 

- -d+d.^+ ^ -0 


(67) 
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where 9± = 2^. 

Defining the Misner-Sharp mass for this metric as M = (1 + [55], the apparent horizon 

for this metric is located in 


R = 2M => d+Rd-R = 0 ( 68 ) 

The future and past apparent horizon are distinguished by = 0 and 9- =0. In the case of gravitational 
collapsing models we consider the future apparent horizon. If we write the general spherically symmetric 
matter tensor in the double null coordinates (u, v) as 

+ T++dudu + _ dudv + T _ dvdv (69) 

then Einstein’s equation give 

Rd+d-R + //2 + d+Rd-R = AnR^T+_. (70) 


Using this result in the general wave equation we get 

f attr^t+_- f/2-d+Rd-R e{e + i)f\ 

2R^) 


We can change the coordinates to {U, V) where metric becomes 

ds^ = -2FdU dV + R^dn? 


(71) 


(72) 


which like the Kruskal coordinates are regular at the horizon and cover all spacetime. The coordinate 
transformation between these two coordinates gives 


f -F(^^ dU dU dV 
^du dv dv du 


{72,) 


Without loss of generality we take U = p{u) and V = q{v). Since the {u,v) coordinates only cover the 
exterior of the horizon and u goes to infinity at the future apparent horizon, we get 

^ /^O. (74) 

An easier way to see this is that one can write the general spherically symmetric metric in the Kodama 
foliation 

ds^ = -e’^dU(l - + y-^^^^ + R{t,r)^dn^ (75) 

R > 

Here M is the Misner-Sharp mass and these coordinates cover outside the horizon. This metric can be 
written 


Then we have 


ds^ = -e’^(‘’’’)(l- 


‘^M{t,r)+^ 


R 


)\dt^ 


dr^ 


( 1 - 


2M{t,r 


-) -I- R{t, r)^dD^ 


/ = /(t,r)e^(‘-’')(l 


2M{t,r) ^ 
R ’ 


dudv = —- vidt^ + — 

(1 


dr'^ 

2M(t,r) ^2-* 

R ) 


(76) 


(77) 


where I{t,r) is an integration factor to ensure that the second term is a perfect differential. One can 
easily see that / is zero on the horizon. Therefore, even in the general spherical symmetric model we can 
neglect the backscattering, £ 7 ^ 0 term in eqn. 0 , and only consider the s-waves term near the future 
apparent horizon, given by 9+ = 2^ = 0. 


Unlike the Schwarzschild black hole where the general wave equation near the apparent horizon reduces 
to the 2 -dimensional wave equation and only the first term in equation 0 takes part, generally non-zero 
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matter terms cause that the 2-dimensional wave equation can be different from the 4-dimensional wave 
equation near the future apparent horizon for generally spherically symmetric models. As a result of the 
above calculation, in the case that we have a zero matter term, we can neglect the effective potential term, 
the 4-dimensional general wave equation reduces to the 2-dimensional wave equation, and the general 
features of the radiation like its thermal character and temperature are the same. On the other hand 
we have discussed above that only isolated black hole (black holes which already have no flux or matter 
on the horizon) can radiate. Therefore, we can infer that the backscattering term can be neglected for 
all radiating black holes, and s-waves and the 2-dimensional wave equation can well describe the general 
features of the radiating black hole near horizon even in the non-static case. 


VI. CONCLUSION AND DISCUSSION 

We have shown that energy conservation, required for virtual pair creation to be real, gives a 
strict limitation on which black holes can be radiative. It was shown that general black holes that 
have a dynamical horizon |40j cannot have pair creation on the horizon, because they violate energy 
conservation. Only isolated horizons and quasi-isolated horizons (slowly evolving horizon) can create 
real pairs which are not forbidden. 

Next, we have considered how the expectation value for the quantum stress tensor carries information 
about quantum effects in curved space time. Beside the zero point energy there are two important 
quantum effects, namely particle creation and vacuum polarization, which can be studied through the 
expectation value of the quantum stress tensor in curved space time. Generally there is no way to 
distinguish these two effects in curved space time, but having a horizon and choosing suitable coordinates 
help us to identify these two effects in the expectation value of the quantum stress tensor. Existence of 
an event horizon so that the observer cannot see part of spacetime is one criterion for having particle 
creation in a space time. To identify particle creation in the quantum stress tensor, one must choose a 
vacuum state that covers both inside and outside the horizon; this is a mixed state for the observer. 

In this way, we have studied particle creation and vacuum polarization in cosmological models. We 
have defined three different vacuum states for the de Sitter spacetime which are like the Boulware, 
Hawking-Hartle, and Unruh vacua in the black hole case. It was discussed that the Boulware like vacuum 
state describes the vacuum polarization in the expectation value of the energy momentum tensor. The 
Hawking-Hartle vacuum state is the state of a thermal equilibrium universe with the temperature 
T = The calculation of the particle creation term by the expectation value of the quantum stress 
tensor gives the same thermal radiation that Gibbons and Hawking have calculated [53] for the de Sitter 
stationary spacetime. It was shown that this thermal temperature is compatible with the temperature 
that is calculated by the affine null parameter approach. In addition, the different components of the 
vacuum polarization tensor are calculated in this frame. 

Different aspects of particle creation in the cosmological de Sitter were discussed. Since we only 
able to see the events which are inside our past light cone, the observable quantities must be written 
for that region. For example in the de Sitter static spacetime we have to calculate the quantities 
for a coordinate that just covers inside the horizon (like cosmological tortoise coordinate). For the 
cosmological de Sitter model, the Unruh like vacuum that we proposed is interesting for cosmological 
application because, like the Bunch-Davies vacuum, it covers also the eternally expanding cosmolog¬ 
ical part (upper triangle) of the de Sitter hyperboloid and gives quantities which an observer can measure. 

For the black hole case, the global nature of vacuum states is discussed and it was shown that it is not 
possible to study the vacuum polarization and particle creation locally by using the expectation value of 
the quantum stress tensor. We have shown that particle flux from the quantum stress tensor gives the 
same temperature that the affine null parameter approach [^ gives for the asymptotically flat spacetime. 
We also investigated the nature of the particle term via Bogoliubov coefficients in the case that the event 
horizon has formed. The exponential relation between affine parameters in future and past null infinity 
implies thermal radiation properties for particles with positive energy and anti-particle with negative 
energy. This shows that the created particles are confined near the event horizon. We have shown that 
the vacuum polarization flux for comoving observer (relative to the t = coordinate) in the case of a 
spherically symmetric static compact star or black hole and de Sitter space time is zero. 
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As an example, using the expectation value for quantum stress tensor, we have shown that the particle 
creation surface cannot be attributed to the IMOTS surface in Oppenheimer-Snyder collapse. Finally, 
it was shown that we can neglect backscattering, £ ^ 0, terms and only consider the s-waves term near 
the future apparent horizon; the s-wave term gives the important part of Hawking radiation in general 
spherically symmetric black holes. 

This paper introduces ways for choosing good coordinates for separating particle creation from vacuum 
polarization in the case of general non-stationary metrics. This will be of practical use in examining the 
black hole back-reaction problem [44] . 

In the Appendices we summarise the nice presentation of quantum stress tensors given by Saida, Harada 
and Maeda m that is used in what is presented in previous sections. Appendix A gives the 2-dimensional 
case and Appendix B the 2-dimensional case. 


Appendix A: Two dimensional case 


It has already been recognized for a few decades that many different methods of renormalization give 
equivalent results (see for example, chapters 6 and 7 in (Hj). We consider a minimally coupled massless 
scalar field (j), whose stress-energy tensor is given by 


T 



(Al) 


The background spacetime is described in double null coordinates [u, v) as 

ds^ =—Diu^v) dudv . (A2) 


The field (p satisfies the Klein-Gordon equation Dp = 0. When a coordinate system (not necessarily 
null) is specified to describe the differential operator □, we can find a complete orthonormal set {fa;} 
for arbitrary solutions of Dp = 0, where w denotes the frequency of the mode function. The positive 
frequency mode is the mode function which is constructed to satisfy the conditions, w > 0 , 
iUJa;') = S{uj - J) , {fa;, ff;>)=0 and if*, f*,) = -6{u} - Uj'), where {f,g) is the inner-product defined 
from the Noether charge of time translation of p and f* is complex conjugate to fa;, called the negative 
frequency mode. 


In two dimensional spacetimes, the positive frequency modes can be decomposed with respect to 
the direction of propagation. In the double null coordinates, they are fa;{u) = exp(—i uj u) / V Aitw and 
faj{v) = exp(—t w'cj/'s/dTra;. Then, the quantum operator p is expanded by the complete orthonormal set 
of the positive and negative frequency modes as 

poo 

P{u,v)= dw [ua; fa;iu) + al fl (u) + ba; fa; (v) + bl /* (v) ] . (A3) 

Jo 

The canonical quantisation presumes the simultaneous commutation relation between p and its conju¬ 
gate momentum, so that {a^j} and {6^^} are harmonic operators satisfying the commutation relations; 
[aa;,al,] = S{uj — oj') and [ba;,bl,] = S{uj — uj') and all others vanish. They define the Fock space of 
quantum states and give particle interpretation. The vacuum state |vac) is defined as a quantum state 
satisfying 0,^1 vac) = 6(^|vac) = 0 for all w. 


If we choose different coordinates {u,v), a natural orthonormal set of mode functions is {fa;}, 
where fa;iu) = exp(—i w u) /V 47rw and fa;{v) = exp{—i ujv)/^/ dTrw. Then the expansion (A31 gives 
different harmonic operators {da;} and {boj}. These harmonic operators define another vacuum state 
|Wc) (^ I vac)) if there arises a mixing of positive and negative frequency modes (/u;,/*/) ^ 0 between 
the two coordinate systems. Thus, even if a quantum state is initially set to be a vacuum state, this does 
not remain vacuum but corresponds to an excited state associated with the coordinate system natural to 
an observer at the final time if the mixing of positive and negative modes arises. This will be interpreted 
as quantum particle creation in curved spacetimes. 
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The quantum expect ation value of the stress-energy tensor (vacjT^y |vac) is calculated by substituting 
the quantum operator (A3) (after replacing and bc^ with and bc^) into the stress-energy tensor 
However, (vac|r^i(|vac) diverges even for flat background cases. Therefore, we need to renormalise 
the stress-energy tensor. We do not go into the details of the regularization method but only quote the 
result §6.4 in [H], 


IT \ -ft -L — 

' AO 3fiy •> 
4o TT 


(A4) 


where (T^p) is the renormalized expectation value of (vac|Tpp|vac), TZ is the Ricci scalar of the background 
spacetime, and is a symmetric tensor whose components in the coordinate system (u, v) on which the 
vacuum |vac) is defined is given by 


h 

'3 

1 

(M 

^ ,uu 

247r 

2 

^ D J 

D 


h 

'3 


D ,vv 

247r 

2 

K D ) 

D 


(A5a) 

(A5b) 

(A5c) 


where D{u,v) = —2guy. The renormalized expectation value (T^p) of (vac|T^i/|vac) in the other coor¬ 
dinates (u, v) is calculated from the above components through the usual coordinate transformation for 
tensor components. 


(T^p) = 


dx^ dx'' 
^x^^ dx'' 




(A6) 


Appendix B: Four dimensional case 


The renormalized vacuum expectation value of the stress-energy tensor {T^^) of matter field in four 
dimensions may also be calculated via the canonical quantization formalism as shown for the two 
dimensional case in the previous section. However the path integral quantization formalism is more 
convenient to summarize {T^^,) on a four dimensional conformal spacetime. 


The effective action IT of a quantum matter field (/) on a spacetime of metric gives the vacuum 
expectation value of a quantum stress-energy tensor. W can be evaluated by the path integral method 
and the vacuum state |vac) is specified by the Green function of (j) used in evaluating the path integral. 
However the precise path integral form of W is not important here. W is decomposed into two parts as 
W = ITren + ITdiv, where Wren is the renormalized part and Wdiv is the divergent part. The functional 
differentiation of Wren gives the renormalized vacuum expectation value (7)rp), 




2 <5Wren 


(HI) 


We consider the case where the metric 5^1/ is conformal to the other one as 

9iJ^y — ^ ^ 9fj-y 5 


(B2) 


and the matter field (() is a conformally coupled massless scalar field satisfying (□ — TZ/6)(j) = 0. On the 
other hand, we get from the definition of functional differentiation. 


Wren “ Wren = 


^Wr< 


Sg°‘P 


5g°‘^ d^x, 


(B3) 


where Wren is the renormalized effective action obtained from Wren on replacing g^y by 'g^y. Then 
considering functional differentiation only by the conformal transformation, Sg^'^ = —2g^'' 5il, the 
effective action is expressed as 

/ xo 

9°"^ {Tap) . (B4) 
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Substituting this into eq. (Bl), we get 






where <5^ = g/j.ag^"', g^ag'^°‘ = 5^ and the general relation, 




6g° 


6g° 


(B5) 


(B6) 


is used to get the first term of the right-hand side of eq. (B5). The trace g^^^ {Tap) is usually called the 
conformal anomaly or the trace anomaly, and it is well known that the divergent part Wdiv gives the 
conformal anomaly as (see §6.3 in m for example) 


g’^^iTap) = 


n ^ITdiv 


6n 


(B7) 


Hence substituting this expression of the conformal anomaly into eq. (B51 and using eq. (B 61 and eq. (B3), 
on replacing Wren by Wdiv we obtain 


(T^u) — ^2 


2 ,5Wdiv 2^2 JWdiv 


y=5 6gf^'' 


(B8) 


The divergent part Wdiv can be evaluated from the Green function of the matter field (p. We do not 
present the details of the calculation of Wdiv, but quote only the result for for the conformally 

coupled massless scalar field ip on the spacetime of metric (B2) (see §6.2 and §6.3 in m for detail). 


{T^v) — ^2 


n 


1 


28807r2 \ 


n 


1 


28807r2f]2 I 6 ' ’ 


where 


■— 2 V^V,y7?. — 2(;^,y □??.-I-- 7?. g^^ — 2TZ , 
- TZ R^i, -I- - RapR°‘^ g^v — 


V := -R^Rcu 




(B9) 

(BlOa) 

(BlOb) 


and R^n and TZ are the Ricci tensor and scalar with resp ect t o g^n respectively, and and are 
defined similarly with respect to the metric Eq uation (B9) is the generalization of eq. (6.141) in [14] 
to the general conformal spacetimes of metric (B 2 ). 
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